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1.  INTRODUCTION 


Curtis,  Srour,  and  Chiu1  recently  studied  the  transport  of  electrons 
and  holes  in  Si02  films  during  exposure  of  the  films  to  ionizing 
radiation.  Electron-hole  pairs  were  excited  in  an  oxide  film  by  a 
pulsed  electron  beam  with  energy  in  the  4 to  8 keV  range.  The 
dependence  of  current  on  applied  electric  field  was  observed  for  fields 
up  to  -5  x 106  V/cm . These  researchers  proposed  that  the  observed 
field  dependence  of  current  could  best  be  explained  by  geminate2  or 
columnar3  recombination  or  both.  By  estimating  the  energy  deposited  by 
the  electron  beam  in  the  Si02  layer  and  assuming  a value  of  27  eV  to 
create  one  electron-hole  pair,  they  concluded  from  their  measurements 
that  substantially  all  of  the  generated  electron-hole  pairs  are 
collected. 

* 

In  the  present  study,  the  data  of  Curtis  et  al1  were  examined  and 
demonstrated  to  be  explained  remarkably  well  by  columnar  recombination. 
A pair-creation  energy  of  about  18  eV  in  Si(>2  has  been  obtained  by 
extrapolation  to  infinite  electric  field  of  the  columnar  model  fit  to 
their  data.  Also,  this  magnitude  of  the  pair-creation  energy  has  been 
independently  verified  by  examination  of  the  energy  loss  mechanisms  of 
low-kiloelectron-volt  electrons  in  Si02 . 

! 

! 

2.  PLASMON  CALCULATION  OF  THE  PAIR-CREATION  ENERGY 

When  a high-energy  electron  traverses  a solid  film  such  as  Si02 , it 
loses  energy  primarily  by  low-momentum  transfer  Interactions  with  the 
valence  electrons.  Plasmon  excitation  is  the  principal  mechanism  of 
low-momentum  transfer  energy  loss  of  fast  electrons  in  solids.1*  The 
plasmon  is  a well-defined  excitation  mode  of  the  valence  electrons  in  a 
solid  that  is  characterized  by  a collective  oscillation  of  all  of  the 
electrons  at  the  "plasma"  frequency.  Although  a portion  of  the  incident 
electron's  energy  may  be  lost  through  the  excitation  of  high-energy 
electrons,  these  electrons  also  lose  their  energy  predominantly  by 
plasmon  production;  hence,  almost  all  of  the  energy  of  the  incident 
' electron  is  ultimately  lost  to  the  plasmon  mode.5  Following  its 

creation,  the  plasmon  decays  rapidly  into  electron-hole  pairs.  For 
determination  of  the  average  energy  to  create  an  electron-hole  pair,  the 
valence  band  structure  of  Si02  must  ba  known.  DiStephano  and 


10.  L.  Curtin,  Jr.,  J.  R.  Srour,  and  K.  Y.  Chiu,  J.  Appl.  Phys.,  45_ 
(1974),  4506. 

2L.  Onsager,  Phys.  Rev.,  54  (1938),  554. 

3G.  Jaffe,  Ann.  Phys.,  42_  (1913),  303. 

I*.P.  Pines , Rev.  Mod.  Phys.,  28  (1956),  184. 

V.  RothwMf,  J.  Appl.  Plxys.,  4£  (1973),  752. 
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Eastman6  determined  the  basic  structure  of  the  Si02  valence  band  by 
photoelectron  spectroscopy.  The  2p  part  of  the  band  consists  of  a 2p 
nonbonding  band  lying  between  0 and  -4  eV  below  the  top  of  the  valence 
band  ,(Ev)  and  a lower- lying  2p  bonding  band  extending  to  - 11  eV  below 
\EV.  The  2p  bonding  band  contains  twice  as  many  electrons  as  the  2p 
nonbonding  band.  The  2s  part  of  the  valence  band  extends  from  -18 
to  -24  eV  below  Ey.  The  schematic  diagram  of  figure  1 represents  the 
approximate  band  structure  of  si02 . 


CONDUCTION  BAND 

Ec~  3«v 


Figure  1.  Schematic  diagram  of  the  Si02  band  structure  adapted  from 
the  work  of  T.  H.  DiStephano  and  D.  E.  Eastman  (Phys.  Rev.  Letters, 
27  (1971) ) . The  approximate  values  delineating  the  various  energy 
bands  are  used  for  schematic  representation  only. 


6T.  H.  DiStephano  and  D.  E.  Eastman,  Phys.  Rev.  Letters,  27_  (1971)  , 
1560. 
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The  plasmon  frequency  is  given  by 


to2  = 4irne2/m  (1) 

P 


where  n is  the  number  of  electrons/cm3  participating  in  the  plasma 
oscillation,  e is  the  electron  charge,  and  m is  the  f ree-electron  mass. 
If  all  the  valence  band  electrons  participated  equally  in  the  plasma 
oscillation,  then  the  plasmon  energy  Ifu  would  be  22.4  eV.  To  produce 
an  electron-hole  pair,  an  electron  in  “the  valence  band  must  be  excited 
into  the  conduction  band  by  receiving  sufficient  energy  to  overcome  the 
band  gap  of  Sit»2  ( eV) . The  decay  of  the  plasmon  yields  sufficient 
energy  to  excite  an  electron  from  any  part  of  the  2p  band,  but  not 
enough  to  excite  an  electron  from  the  2s  band.  When  the  plasmon  decays 
by  exciting  an  electron  from  the  2p  nonbonding  band,  the  excited 
electron  has  sufficient  kinetic  energy  to  excite  an  additional  electron 
from  the  2p  nonbonding  band.  Although  the  average  number  of  pairs 
created  probably  lies  somewhere  between  one  and  two,5  for  the  calcula- 
tions below,  the  plasmon  is  assumed  to  yield  two  electron-hole  pairs 
when  it  decays  by  exciting  electrons  from  the  2p  nonbonding  band.  Since 
the  number  of  states  in  the  2p  bonding  band  is  approximately  twice  that 
in  the  2p  nonbonding  band,  twice  as  many  plasmons  are  expected  to  decay 
by  exciting  an  electron  from  the  2p  bonding  band  as  from  the  2p 
nonbonding  band.  Therefore,  on  the  average,  for  every  three  plasmons 
created,  two  decay  by  exciting  one  electron  each  from  the  bonding  band, 
and  one  decays  by  exciting  two  electrons  from  the  nonbonding  band. 
Thus,  three  plasmons  decay  into  four  electron-hole  pairs.  From  the 
value  22.4  eV  for  the  plasmon  energy,  the  average  energy  required  to 
produce  one  electron-hole  pair  is  estimated  to  be  W * 
3 x 22.4/4  = 16.8  eV.  This  estimate  of  the  pair-creation  energy 
compares  very  well  with  the  value  obtained  below  when  the  columnar  model 
of  recombination  is  applied  to  the  data  of  Curtis  et  al.1 


10.  L.  Curtis  Jr.,  J.  R.  Srour,  and  K.  Y.  Chiu,  J.  Appl.  Phys.,  45 

(1974),  4506. 

5A.  Rothwarf,  J.  Appl.  Phys.,  44  (1973),  752. 
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The  mean  free  path  for  the  plasmon  production  is  given  by  7»8»* 


X = 2 
P 


( ( E/fioip  jy/  In  |4E/fia)p  | 


(2) 


where  a^  is  the  Bohr  radius  (0.53  A)  and  E is  the  energy  of  the  fast 
electron.  By  applying  the  work  of  Everhart  and  Hoff  9 to  the  Al-SiCfc-Si 
system  with  an  Si02  layer  of  1050  A,  Curtis  et  al1  obtained  depth-dose 
curves  for  various  combinations  of  initial  electron  energy  and  SiOg  film 
thickness.  From  these  curves,  am  electron  having  an  initial  energy  of 
4 keV  is  estimated  to  have  an  average  energy  of  approximately  1.5  keV  in 
the  Si02  layer.  For  1.5  keV  electrons,  X = 12.7  X.  When  three 
plasmons  decay  by  creating  four  electron-hole  pairs,  the  mean  free  path 
for  electron-hole  excitation  is  X=  9.52  A.  This  value  for  X,  coupled 
with  the  energy  per  electron-hole  pair,  is  consistent  with  the  energy 
loss  of  an  electron  traversing  the  insulator  thickness  as  determined 
from  the  depth  dose  curves  of  Curtis  et  al.1  These  curves  indicate  an 
energy  loss  of  ~1.8  keV  for  an  electron  traversing  1050  a of  Si02  film. 
On  the  other  hand,  from  the  calculated  values  of  X and  W,  the  energy 
loss  across  1050  A of  Si02  is  estimated  to  be  (1050  A/X)W  = 1.85  keV. 

The  small  value  of  X,  the  separation  between  adjacent  electron-hole 
pairs,  justifies  the  columnar  model  of  ionization. 


^0.  L.  Curtis,  Jr.,  J.  R.  Szouz , and  K.  Y.  Chiu,  J.  Appl.  Phys.,  45_ 
(1974),  4506. 

7 D.  Pines,  Phys.  Rev.,  92_  (1953),  626. 

8 A.  Rose,  RCA  Review,  27_,  No.  4 (December  1966),  600. 

9 T.  E.  Everhart  and  P.  H.  Hoff,  J.  Appl.  Phys.,  42_  (1971),  5837. 

* Actually,  this  expression  for  Xp  includes  energy  loss  to 

single-electron  excitations , as  well  as  plasmon  excitations.  Therefore, 
equation  (2)  underestimates  the  mean  free  path  for  plasmon  excitation  by 
the  incident  high-energy  electron.  However,  the  secondary  electron 
produced  in  a single-electron  excitation  loses  energy  in  the  same  way 
as  the  incident  electron.  Following  a single-electron  interaction,  the 
secondary  electron  and  the  incident  electron  both  have  less  energy  than 
the  initial  energy  of  the  incident  electron.  Also,  Xp  is  approximately 
proportional  to  the  electron  energy.  Therefore , equation  (2)  still 

represents  a reasonably  good  estimate  for  the  effective  mean  free  path. 
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3.  THEORY  OF  COLUMNAR  RECOMBINATION 


The  columnar  (or  track)  model,  developed  by  Jaffe,3  was  originally 
applied  to  ionization  in  gases  at  high  pressures.  The  columnar  theory 
is  based  on  a cylindrically  symmetrical  distribution  of  ions  about  the 
path  of  the  ionizing  particle  as  axis , the  density  of  ionization  being 
at  all  points  the  same  for  positive  and  negative  ions  and  falling  off 
with  increasing  distance  from  the  axis.  The  ordinary  laws  of  diffusion, 
bimolecular  recombination,  and  drift  under  an  electric  field  are  applied 
to  this  distribution,  and  the  fraction  of  ionized  pairs  escaping 
recombination  is  calculated.  The  application  of  the  macroscopic  laws  of 
bimolecular  recombination  tacitly  assumes  that  the  mean  separation  of 
charges  of  opposite  sign  exceeds  the  critical  radius  for  recombination. 
(The  critical  radius  is  that  distance  of  separation  of  a positive  and 
negative  ion,  within  which  recombination  is  certain  and  outside  which 
escape  is  probable.)  Also,  the  neglect  of  geminate  recombination?  (the 
recombination  of  an  electron  with  the  same  molecule  from  which  it  was 
ionized)  is  based  on  the  assumption  that  ionization  along  the  path  of 
the  primary  ionizing  particle  is  so  dense  that  there  is  no  correlation 
of  ionized  electrons  with  their  parent  ion  and  that,  therefore, 
homogeneous  bimolecular  recombination  occurs  within  the  column. 

O 

In  his  model  for  columnar  recombination,  Jaffe  assumed  that  both 
the  positive  and  negative  charge  species  have  the  same  mobility.  In  the 
present  study,  his  model  was  modified  to  allow  the  two  charge  species  to 
have  different  mobilities.  This  modified  model  should  be  more 
appropriate  for  a solid  such  as  Si02  , in  which  the  mobility  of  holes  is 
probably  much  less  than  that  of  electrons.  The  current-field 
relationship  resulting  from  the  modified  Jaffe  theory  is  given  by 


where 


1 = 


1 + joN  /4ir(Di  + D2>!  e*  K (x) 
| o ) o 


x = (y i + (j2)2<sin2  0>  b2  E2/  j4  (Dj 


(3) 


(4) 


and  K is  the  modified  Bessel  function  of  order  zero.  In  this 
relationship,  a is  the  recombination  coefficient,  N is  the  linear 
ionization  density  (N  = A-1),  Di  and  D2  are  the  diffusion  coefficients 
for  electrons  and  holes,  Uj  and  u2  are  the  mobilities  of  the  electrons 
and  holes,  <sin2  0>  is  the  average  over  a random  distribution  of  the 


2L.  Onsager,  Phys.  Rev.,  54.  (1938),  554. 

3G.  Jaffe,  Ann . Phys.,  42.  (1913),  303. 
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angle  between  the  electric  field  E and  the  column  axis,  b is  a measure 
of  the  column  radius,  and  I is  the  saturation  current.  At  high  fields. 


I/IQ  » l/£l  + (a  NQ/yb)  >fi72/4E^  . 


(5) 


Thus,  at  sufficiently  high  fields,  a plot  of  I versus  E should  yield 
a straight  line  with  intercept  I -1  . Since  a and  y are  related  by  the 
Langevin  relation  1 0 


a = |e/eQejy  (6' 

(where  y = y j + y ^ , and  N = X “ 1 has  been  estimated  from  the  plasmon 
calculation,  the  slope  of  °the  straight-line  plot  of  I-1  versus  E-1  at 
high  fields  gives  a direct  measure  of  b,  che  column  radius.  The  direct 
dependence  of  equation  (3)  on  a,  D,  and  y is  eliminated  by  the 
Langevin  relation,  equation  (6) , and  the  Einstein  relation  between  y and 
D.  The  only  free  parameters  of  the  theory  are  ID,  b,  and  X.  I0  and  b 
are  determined  from  the  least-squares  fit  to  the  high-field  data, 
and  A is  estimated  from  the  plasmon'-energy-loss  calculation. 


4.  DISCUSSION  OF  RESULTS 

° i i 

From  the  data  on  the  1050  A sample  of  Curtis  et  al , * I 1 versus 
E ~'has  been  plotted,  and  good  straight  lines  have  been  obtained  in  the 
high  field  region.  From  the  intercept  (I  -1),  the  values  of  electron 
beam  current  ( 1^)  and  beam  energy  (Eg),  and  the  estimated  fraction  of 
energy  deposited  in  the  SiO  2,  the  pair-creation  energy  w can  be 
obtained.  For  E0  = 4 keV  and  I - 0.52  nA,  a least-squares  fit  to  the 

high  field  data  yields  I = 60.5  nA,  b = 65.6  &,  and  W = 16.4  eV.  This 
value  of  W compares  very  well  with  our  estimate  based  on  plasmon 
generation  and  subsequent  decay  into  electron-hole  pairs.  With  these 
values  for  I and  b and  the  value  A = 9.52  A (N  = A-1) , the  complete 
field  dependence  of  I can  be  calculated  from  equation  (3) . Figure  2 
shows  the  excellent  agreement  between  the  experimental  values  of  current 
and  the  values  calculated  from  the  columnar  model. 


10.  L.  Cuztis,  Jr.,  J.  R.  Sronr,  and  K.  Y.  Chiu,  J.  Appl.  Phys.,  45_ 
(1974),  4506. 

13P.  Langevin , Ann.  Chim.  Phys.,  28_  (1903),  289,  433. 
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APPLIED  FIELD  (V/cm) 

Figure  2 . Voltage  dependence  of  electron-beam  irradiatic.i-induced 
current  in  Si02.  The  circles  (from  0.  L.  Curtis,  Jr.,  J.  R.  Srour, 
and  K.  Y.  Chiu,  J.  Appl.  Phys.,  £5  (1974))  represent  data  from  a 
1050  A sample  irradiated  by  a 4 kV  electron  beam  with  a beam  current 
of  0.52  nA.  The  solid  line  is  the  result  of  the  columnar  model 
calculation. 


For  values  of  W and  b averaged  over  different  sets  of  experimental 
data,  the  Jaffe  model  was  applied  to  a normalized  current  versus  field 
curve  representing  data  for  several  different  beam  currents.  The  values 
W = 18.4  eV  and  b = 81,3  A were  obtained  from  this  analysis. 

After  the  analysis  described  in  this  report  was  carried  out,  Srour 
et  al1 11  reported  additional  experiments  on  Si02  films,  which  were 
prepared  by  a different  processing  technique.  Carrier  collection  at 
high  electric  fields  increased,  presumably,  due  to  a reduced 
concentration  of  recombination  centers  in  these  oxides  grown  by  the 
different  process.  Srour  et  al11  concluded  that  the  carrier  collection 
at  the  higher  fields  establishes  an  upper  limit  of  ~ 19  eV  for  the 
electron-hole  pair-creation  energy. 


1 R.  Srour,  0.  L.  Curtis,  Jr.,  and  K.  Y.  Chiu,  IEEE  Trans.  Nucl. 

Sci.,  NS-21  (December  1974),  73. 
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Earlier  in  this  report,  it  was  assumed  that  plasmon  decay  results  in 
two  electron-hole  pairs  excited  from  the  2p  nonbonding  band.  Since  the 
average  number  of  pairs  created  from  this  band  per  plasmon  probably  lies 
somewhere  between  one  and  two,  the  value  W = 16.8  eV,  obtained  from  the 
calculation  of  plasmon  generation  and  subsequent  decay  into  electron-hole 
pairs,  represents  a lower  limit  on  the  pair-creation  energy.  The  upper 
limit  on  W is  the  full  plasmon  energy,  22.4  eV.  Similar  arguments  require 
that  the  value  X = 9 . 5 A be  a lower  limit  for  the  mean  free  path  for 
electron-hole  excitation.  The  actual  value  of  X probably  lies  between 
9.5  and  12.7  A,  the  mean  free  path  for  plasmon  production.  The  only 
important  effect  on  the  columnar  model  of  a different  value  for  X is  to 
change  the  value  of  b through  the  relationship  Xb  = constant.  The 
calculated  values  of  1 should  be  unaffected,  except  possibly  at  the 
lowest  fields,  since  at  moderate  to  high  fields  the  columnar  model 
depends  only  on  the  product  Xb. 
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